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ABSTRACT 

In the context of recent observational results that show massive ellipticals were in place at high 
redshifts, we reassess the status of monolithic collapse in a ACDM universe. Using a sample of over 
2000 galaxies from the Sloan Digital Sky Survey, by comparing the dynamical mass and stellar mass 
(estimated from colours) we find that ellipticals have ‘cores’ which are baryon-dominated within their 
half-light radius. These galaxies correspond to 3-sigma peaks in the spherical collapse model if the 
total mass in the halo is assumed to be 20 times the dynamical mass within the half-light radius. This 
value yields stellar mass to total mass ratios of 8%, compared to a cosmological baryon fraction of 
18% derived from WMAP3 alone. We further develop a method for reconstructing the concentration 
halo parameter c of the progenitors of these galaxies by utilizing adiabatic contraction. Although 
the analysis is done within the framework of monolithic collapse, the resulting distribution of c is 
log-normal with a peak value of c ~ 3 — 10 and a distribution width similar to the results of N-body 
simulations. We also derive scaling relations between stellar and dynamical mass and the velocity 
dispersion, and find that these are sufficient to recover the tilt of the fundamental plane. 

Subject headings: dark matter - galaxies: elliptical and lenticular, cD - galaxies : fundamental pa¬ 
rameters - galaxies: evolution - galaxies: formation 


1. INTRODUCTION 

The la rge new galaxy surveys such as 2dF 
JColless et, all 120011) and the Sloan Digital Sky Survey 
(SDSS) itYork et a,1.1120001) have transformed the way we 
can study galaxy properties statistically. Here we focus 
on a sample of elliptical galaxies derived from those 
selected from the SDSS bv lBernardi et ali l|2003al) (here¬ 
after B03). Our motivation is to revisit fundamental 
issues of galaxy formation from the perspective afforded 
by a modern data set. 

A theory of galaxy formation in a universe dominated 
by Cold Dark Matter was first set out in a seminal paper 
by Blumenthal, Faber, Primack & Rees (1986a) (here¬ 
inafter BFPR) which considers the monolithic collapse 
of isolated dark matter halos. According to this picture, 
the redshift at which a halo collapses is determined only 
by its mass, for a given cosmology and choice of ampli¬ 
tude of fluctuation. The baryons follow the dark mat¬ 
ter distribution until the radius of the collapsing halo 
reaches the virial radius and is halted. Studies of the 
evolution of the fluctuations in density, however, suggest 
that galaxy formati on is do minate d by merging of small 
dark matter halos iWhite & Rees LL978D, and this pic - 
ture is supported by simulations I Snringel et al • 2QflJ5). 
There have been many attempts to discriminate between 
these two pictures of the formation of galaxies, referred 
to respectively as monolithic or spherical collapse and hi¬ 
erarchical merging. Simulation results notwithstanding, 
the small scatter of the observed colour magnitude re¬ 
lation and its evolution with redshift provides evidence 
that massive ellipticals were already in place at a red¬ 
shift of z ~ 1 — 2 with little subsequent merging (e.g. 
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iStanfonLetjaD fHHf)), whereas less massive ellipticals 
prese nt features characterist ic of recent star formation 
(e.g. iFerreras fc Silkl (| 2000aD). This ‘inve rted hierarchy’ 
or ‘downsizing’ effect llCowie et al. 1 119961) illustrates the 
complexity of galaxy formation compared with a simple 
picture of the assembly of dark matter halos. In this 
paper we use the predictions of the spherical collapse 
scenario as a benchmark, which can then be challenged 
with further comparisons with observations and detailed 
simulations. 

A new aspect of our analysis is the estimation of the 
stellar mass of each of the B03 ellipticals, allowing us to 
distinguish between baryonic and total mass. We also 
update the calculations of BFPR for ‘concordance cos¬ 
mology’ (matter density f l m = 0.3,dark energy density 
Ha = 0.7,Hubble constant h = 0.7 and amplitude of 
fluctuations = 0.9) 3 and place the observed ellipti¬ 
cals on a revised cooling diagram. Despite the attention 
BFPR has received, this relatively simple generalization 
to ACDM has not previously appeared in the literature. 
Our analysis of the B03 galaxies also appears to be the 
first attempt to incorporate data for individual galax¬ 
ies (as opposed to schematic data) in this parameter 
space. From the first three years of WMAP observations 
of the cosmic microwave background alone (Spergel et 
al. 2006), H m = 0.24to'o4 and fl b = 0.042±°;°° 3 . For 
simplicity, we use a ‘cosmological’ baryon fraction of 1/6 
unless otherwise sta ted. 

On galactic scales, IKlvninl ll200l model the Milky Way 
within the virial radius (« 250 kpc) and find that a 
substantial ‘feedback’ mechanism which removes baryons 
from the galaxy itself must have operated. Do similar 
processes operate in massive elliptical s? Are their cores 
domin ated by baryonic matt er? iRomanowskv et aU 
120031) and iDekel et~ahl ll2005f) have reached conflicting 

3 We note that the recent WMAP release prefers Dm ~ 0.24 and 
<78 ~ 0.74. However, <78 = 0.9 is still preferred by weak lensing. 
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conclusions to this important question, based on the 
analysis of planetary nebulas in elliptical galaxies out to 
five times the effective radius R e . By using the central 
velocity dispersion and stellar mass out to R e we seek 
to answer the same question by focussing on the central 
region only. We find that the cores of the galaxies are 
baryon dominated, at least to a distance of R e from the 
centre. 

By making the straightforward assumption that the to¬ 
tal mass of the galaxy is proportional to the mass of the 
galaxy ‘cores’ we have studied, we reproduce the BFPR 
results for a modern ‘concordance’ cosmology, and in¬ 
vestigate the regime in which their benchmark model 
of spherical collapse is consistent with the data. We 
also recover the profile of the undisturbed dark matter 
halo from present-day observables via the procedure of 
adiabatic contraction, which enables us to compare our 
derived profile concentration with both simulations and 
observations. Defining a = M tot / Md yn (< R e ) - the to¬ 
tal mass of the collapsing halo divided by the dynamical 
mass within R e - we find a strong constraint that a > 10. 
a « 20 provides a good fit to the data. 

The outline of the paper is as follows: Section 2 
presents the B03 sample and a volume-limited subset of 
2040 galaxies selected from it. Section 3 describes the 
derivation of stellar masses of these ellipticals from the 
SDSS colours. Section 4 derives dynamical masses and 
Section 5 covers the relations between luminosity, ve¬ 
locity dispersion, s tellar mass and dynami cal mass, ex¬ 
tending the classic iFaber fc Jackson ( 197 61 relation. In 
Section 6 we comment on the acceleration of stars in 
the ellipticals and the connection to MOND. Section 7 
extends the benchmark spherical collapse model for the 
concordance A-CDM cosmology and use it to contrast 
the observed ellipticals with the level of rms fluctuations 
on the cooling diagram. Finally, in Section 8 we present 
a new reconstruction method which recovers the initial 
halo of dark matter from the present-day observables us¬ 
ing the adiabatic contraction. We summarize the results 
and outline future work in Section 9. 

2. THE SAMPLE 

We use ga l axies selected from the sample of 
iBernardi et al.! (l2003afl extracted from the SDSS. 4 The 
sample comprises almost 9,000 early-type galaxies, with 
spectroscopic redshifts in the range 0.01 < z < 0.3 (in 
the heliocentric reference frame). In B03, the authors 
describe the selection criteria for their sample and the 
catalogue itself; as the photometry in each of the bands 
is used in our subsequent analysis, it is important to note 

4 Funding for the creation and distribution of the SDSS Archive 
has been provided by the Alfred P. Sloan Foundation, the Partic¬ 
ipating Institutions, the National Aeronautics and Space Admin¬ 
istration, the National Science Foundation, the U.S. Department 
of Energy, the Japanese Monbukagakusho, and the Max Planck 
Society. The SDSS Web site is http://www.sdss.org/ The SDSS 
is managed by the Astrophysicai Kesearch Consortium (ARC) for 
the Participating Institutions. The Participating Institutions are 
The University of Chicago, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, The Johns Hopkins Univer¬ 
sity, the Korean Scientist Group, Los Alamos National Laboratory, 
the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck- 
Institute for Astrophysics (MPA), New Mexico State University 
University of Pittsburgh, University of Portsmouth, Princeton Uni¬ 
versity, the United States Naval Observatory, and the University 
of Washington. 



Fig. 1.— Stellar mass function (MF) of our volume-limited sam¬ 
ple of early-type galaxies extracted from Bernardi et al. (2003). 
As shown in the insert (in which we plot the full sample of ~ 8700 
galaxies), we initially make cuts in redshift and absolute magni¬ 
tude, plotting the MF for the remaining galaxies as the jagged 
line. In order to check that these cuts produce a useful sample, 
we compare the result with the luminosity function derived for 
all SDSS early-type systems (type 0 < T < 1) obtained indepen¬ 
dently by Nakamura et al (2003). For the purposes of this compar- 
ision, we assume a characteristic M s tar/Lr = 4 (see figure Our 
sample is complete for systems with stellar masses greater than 
7.94 X 10 10 Mq (dashed line), and so we make a final cut at this 
mass to leave our final sample of 2040 galaxies. 


that these criteria do not include selection by colour. 
Their sample is magnitude-limited with all galaxies at 
a redshift of z > 0.3 excluded. For the work presented 
in this paper we have used the data as presented in B03, 
with the exception of photometry where we have used the 
recen t Sloan Data Release 4, iAdelman-McCarthv et alJ 
12005 ) which resolves many of the problems with photom¬ 
etry identified in earlier releases. 

Unlike the work presented in B03 and its companion 
papers, we have further restricted the sample to provide a 
volume-limited sample. The mass function of the subset 
of galaxies selected by B03, with z < 0.1 and absolute 
magnitude in the r-band M r < —20.75, is shown in figure 

IH 

As the figure shows, restricting our sample to massive 
elipticals, with M star (< R e ) > 7.94xl0 10 M Q produces 
a volume limited sample of 2040 galaxies with z < 0.1 
which appears free from selection effects. Unless oth¬ 
erwise stated, this is the sample we use for subsequent 
results. The definition of stellar mass, as shown in the 
figure, is given in the next section. 

3. STELLAR MASS 

For massive elliptical galaxies with little gas content, 
the baryonic mass is approximately equal to the stel¬ 
lar mass at least within the effective radius, R e . We 
compute the stellar mass content within the effective ra¬ 
dius by comparing the available SDSS photometry with 
a wide range of star formation histories (SFHs). We as¬ 
sume an exponentially decaying star formation rate, so 
that each SFH is uniquely defined by three parameters: 
stellar metallicity (constant with time), star formation 
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timescale and formation redshift (allowed in the range 
10 > zf > 2). These models can represent either a 
short, sharp burst of star formation or a more extended 
episode. For every choice of SFH we cal culate the result¬ 
ing synthetic stellar populations of iBrnzual fcCliariot 
(2003), generating a composite spectrum. A iChabrien 
(2003) initial mass function is used throughout. We use 
four SDSS passbands to compare the flux in each of three 
colours (g — r,r — i,i — z) with the observations of each 
source. The best fit is used to determine the stellar mass- 
to-light ratio which allows us to estimate the stellar mass 
content of each source from the absolute r-band magni¬ 
tudes given in data release 4. 

It seems initially that the use of the full spectral energy 
distribution (SED) would produce a stellar mass content 
which might be more accurate. However, we must em¬ 
phasize that our results indicate that the stellar mass to 
light ratio in the r band does not change much with re¬ 
spect to the star formation histories allowed in early-type 
galaxies. All photo-spectroscopic observables of local 
early-type galaxies correspo nd to o ld, pas sively evolving 
stellar populations llBernardi et al.ll2003bll . which results 
in a very weak dependence of the M/L for the allowed 
ages and metallicities (see the right hand panel of figure 
0. The largest systematic error in the stellar mass is due 
to the choice of the initial mass function (IMF), which 
amounts to a factor ~ 0.1 — 0.2 d ex in l og (M/L ) for 
a rang e of standard IMFs (see e.g. iBruzual fc Charlotl 
(»)). 

4. DYNAMICAL MASS 

We can use the projected size and central velocity dis¬ 
persion along with a simple model of the dynamics of 
early-type galaxies to estimate the total matter content. 
Assuming that the velocity dispersion and the anisotropy 
of the velocity distribution do not vary w ith radius, one 
can solve .leans’ equation (see e.g. iBinney & Tremaine 

cm)). 

Padmanabhan et al. (2004) perform a careful analysis 
of the profiles of massive galaxies and find that 

M dyn {<R e )=A (1) 

where R e is the effective (or projected half-light) radius, 
a is the velocity dispersion as measured at the centre and 
A is a constant which depends weakly on the anisotropy 
of the velocity field. We adopt the Padmanabhan value of 
A ss 2.72 throughout the present work, noting that they 
assume a systematic error of 30%. The dominant source 
of scatter is, as Padmanabhan et al. note, measurement 
errors. 

Figure [5] presents a comparision of the dynamical mass 
and stellar mass. It is clear that dark matter is under- 
abundant (compared to the cosmological abundance) in 
the central regions of these systems. In other words, the 
core region is baryon dominated, in agreement with ob¬ 
servations (e.g. Mamon & Lokas 2005). The slope is in 
agreement with independent estimates of total mass us¬ 
ing strong gravitational lensing (Ferreras et al. 2005) and 
suggests that dark matter is more important in massive 
galaxies. It also is comparable to the results of Loewen- 
stein and White (1999), although they use an unrealistic 
model for anisotropy in their sources. 



Fig. 2.— Ratio of dynamical to stellar mass measured within R e . 
The solid and dashed lines give the mean and RMS of our 2 : < 0.1 
sample. The dotted line labelled CMB+LSS gives the prediction 
of the cosmological ratio. The inset shows a histogram of the first 
and fourth quartiles in stellar mass (i.e. log (M s tar/M q < 10.72 
and > 11.02, respectively). 


Romanowsky et al. (2003) measure a through observa¬ 
tions of bright planetary nebulae in three local ellipticals 
and derive a mass to light ratio within 5 R e compatible 
with those of purely stellar populations. Dekel et al. 
(2005) explain this result by suggesting that ‘the stellar 
orbits in the outer regions ... are very elongated’. Our 
results, derived for the innermost R e suggest that results 
similar to those obtained by Romanowsky et al. would 
be expected even in large samples of ellipticals. In addi¬ 
tion, realistic assumptions about the velocity strucutre of 
the cores of these systems do not provide an explanation 
of the deviation observed from the cosmological baryon 
fraction. 

5. SCALING RELATIONS 

Figure [3] shows the scaling relations of the (r-band) 
luminosity, stellar and dynamical mass as a function of 
velocity dispersion. The bott om panel shows the we ll- 
known Faber-Jackson relation llFaber fc Jacksoiilll976l) . 

Notice that the truncation imposed by the choice of 
a volume-limited sample results in a constraint nearly 
aligned with the correlation (dashed line). Hence, a 
strong bias is expected if we use this reduced sample for 
the fits (see e.g. appendix B in Lynden-Bell et al. 1988). 
Therefore we use the complete sample for the fits in this 
figure. 

In order to reduce the effect of outliers in the calcula¬ 
tion of the best fitting parameters, in common with all 
fits in this paper we use a robust M-estimator based on 
the minimisation of the mean absolute deviation (Press 
et al. 1992), as standard least squares methods were 
found to be very sensitive to outliers in this sample. 
Given the uncertainties in determining the errors in ve¬ 
locity dispersion and radius, we have chosen to weight all 
points equally rather than take into account their quoted 
errors. A bootstrap method is used to determine the un¬ 
certainty in the slope. The same fitting technique was 
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Fig. 3.— Correlations between velocity dispersion and luminos¬ 
ity (bottom), stellar mass ( middle ) and dynamical mass (top), all 
measured within R e . Characteristic error bars are shown in each 
panel. The solid lines are fits to the data, with all points equally 
weighted. The truncation of our volume-limited sample is shown 
as a dashed line. Given that the truncation and the fit are nearly 
aligned, we have used the complete sample for the fits shown in 
this figure. 


Fig. 4.— Stellar (top) and dynamical mass-to-light ratios as 
a function of velocity dispersion for our volume-limited sample, 
after cuts have been applied. The straight lines are fits to the 
data, and the corresponding power law indices are labelled in each 
panel. The panel on the right shows the predictions for simple 
stellar populations from the models of Bruzual Sz Chariot (2003) 
as a function of age. The dashed, solid and dotted lines correspond 
to metallicities 2, 1 and 1/2 times solar. The histogram refers to 
the data (left panel), with the average given as a horizontal line. 


applied on 1000 bootstrapped samples of 500 galaxies 
each. The predicted slopes from each run were used to 
determine the quoted uncertainties. 

The Faber-Jackson relation is a direct consequence 
of the dynamics of the systems being studied. Early- 
type galaxies are hot dynamical systems whose support 
comes mostly from the velocity dispersion of the con¬ 
stituent stars. Hence, a correlation is expected between 
the central velocity dispersion and the mass of the sys¬ 
tem. The different slope between stellar and dynamical 
mass estimates is mostly related to either stellar pop¬ 
ulation effects, to structural differences, or to a signifi¬ 
cant correlation between the dark matter fractional con¬ 
tribution and total galaxy mass. This effect has been 
commonly invoked to explain the tilt of the observed so- 
called ‘fundamental plane’ with r espect to the predictio n 
from the virial t heorem (see e.g. iGuz nmn^eta D jHH; 
iFerreras fc SilU (l2000bf) : iTniiillo et all 112004) 1. 

Figure^Jshows the correlation between velocity disper¬ 
sion and M/L (stellar or total). Local early-type galax¬ 
ies are dominated by stars which are older than about 
7 — 8 Gyr, which implies that the contribution of stel¬ 
lar populations to the tilt of the fundame ntal p l ane can¬ 
not b e too large, as previously suggested llDressler et all 

GUI- 

There is a remarkable difference between stellar and 
dynamical M/L in figure QJ The fits were done using a 
bootstrap minimisation of the mean absolute deviation 
on the volume-limited sample. The slope of the dynami¬ 
cal M/L is sufficient to explain the tilt of the fundamen¬ 
tal plane. Notice that equation Q] assumes homology for 
all galaxies. Hence, we recover the tilt without invoking 
structural non-homology. Furtherm ore, recent estimates 
using strong gravitational lensing llFerreras et 
suggest dark matter is the main cause of the tilt, a view 


supported by these results. 

6. ACCELERATION 

Our analysis so far has been done within the frame¬ 
work of Newtonian dynamics. We showed that the cores 
(i.e. within radius R e ) of massive early-type ellipticals 
are dominated by by baryons, rather than by dark mat¬ 
ter. The alternative model o f Mod ified Newtonian Dy¬ 
namics (MOND) of lMilgrom lll983i) suggests a modifica¬ 
tion to the acceleration if it is below ao = 10~ 8 cms -2 . 
This modification can then explain the rotation curves 
of spiral galaxies without the need for dark matter. It is 
intriguing indeed to see a correlation between the New¬ 
tonian M/L and the centrifugal acceleration (o = v 2 /r ) 
at the last meas ured radia l point r in spiral galaxies 
( Sanders fc McGaugh 12002^ . Their points are repro¬ 
duced in our Figure 5 (as circles). It is clear from the 
start that such low acceleration is not expected at the 
cores of ellipticals which have much smaller radii than 
the outer parts of spiral galaxies. However, it is interest¬ 
ing to look for any such correlation in the sample of el¬ 
lipticals. Figure 5 shows (in dots) the dynamical M/Lk 
vs. the acceleration, defined as a = v 2 /R e = 3cr 2 /R e . 
We see no obvious correlation and no continuity between 
the spirals and the ellipticals. As expected the ellipticals 
have accelerations which are ~ 10 times larger than that 
required by MOND to show departure form the Newto¬ 
nian dynamics. If indeed ellipticals have large halos it 
would be an important test to measure the acceleration 
at these large radii. 

7. SPHERICAL COLLAPSE IN A ACDM UNIVERSE 

We have already mentioned the renewed interest in 
spherical collapse models as part of the complexity of 
galaxy formation. How do the predictions of this sim- 
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Fig. 5. — Comparison of the acceleration of early-type and late- 
type galaxies. We plot the dynamical M/Lk as a function of accel¬ 
eration (defined as a = u 2 / R e = 3 cr 2 /R e . In order to generate K 
band M/L ratios, we use r — K = 2.4, typical of early-type systems. 
The late-type data (circles) come from Sanders Sz McGaugh(2002). 
Notice that early-type galaxies fall outside of the low-acceleration 
regime required by MOND to give a significant departure from 
Newtonian dynamics. 


Fig. 6.— A comparision of the spherical collapse model with 
our sample of massive early-type galaxies for both EdS (dashed 
lines) and ACDM cosmologies. (EdS parameters : = 1.0,£2 a = 

0.0 ,h = 0.7,cr8 = 0.9, ACDM parameters : Qrn = 0.3,Oa = 0.7,h = 
0.7, erg — 0.9). In obtaining Mtot we have assumed a = 20 (see 
text for details). The Mtot obtained are consistent with a formation 
process from 2- or 3-sigma perturbations in a ACDM cosmology. 


pie model compare to the modern data we have avail¬ 
able, and how do the predictions change in a cosmologi¬ 
cal model incorporating a cosmological constant? There 
have been many attempts to study spherical collapse in 
ACDM (e.g. Lahav et al. 1991; Wang & Steinhardt 
1998) but we will follow the analysis of BFPR more 
closely, using the velocity dispersion-mass and baryon 
density-temperature plots to compare the predictions of 
the model with data. 

In order to check whether a simple spherical collapse 
model is compatible with the properties of massive galax¬ 
ies, we show in figure ED the correlation between total 
mass (i.e. including the halo) and the central velocity 
dispersion. The lines represent the theoretical mass and 
velocity dispersion of dark matter haloes produced by 
monolithic collapse. The collapse of a spherical, top-hat, 
perturbation is described by: 


where D (z) is the growth factor from linear theory. 

Collapse to a single point in this manner is prevented 
by violent relaxation (‘pha se mixing’) which produc es a 
state of virial equilibrium. iBrvan fc Norma.iil ( 1998 1 use 
the virial theorem to obtain the final over-density (rel¬ 
ative to the critical density at the collapse redshift) for 
universes in which fl m + fl A = 1: 

A c = 18 tt 2 + 82 d - 39 d 2 (4) 

where d = fit — 1 and 


m n m {i + zf+ n A + n k (i + z) 2 ' 

A halo will therefore collapse at a redshift, z, determined 
only by its total mass and by cosmolo gy. Via the virial 
theor em its radius will be given by IBarkana fc Loebl 
12001 : 


d 2 r 

df 2 


Hfo A r 


GM (< r) 
r 2 ’ 


( 2 ) 


where r is the radius in a fixed coordinate frame, M (< r) 
is the mass enclosed within r and the initial velocity field 
is equivalent to the Hubble flow only. If the mass de¬ 
scribed by such an equation is bound, it will expand ini¬ 
tially before reaching a maximum radius and collapsing. 
Linear theory predicts that the over-density at whic h the 
top hat collapses to a point is (Wear = 1.686 ( Peebles! 
ft . This value is appropriate for an Einstein-de Sit¬ 
ter universe, but the value is only slightly sensitive to 
the choice of cosmology. Thus, for a top hat to be on 
the point of collapsing at a redshift z, its density extrap¬ 
olated to the present day must satisfy: 


S c (z) 


1.686 

dV) 


= 0.784 


M tl 


10 8 h~ 1 M Q/ 


y 


) 

Tm 1871-2 - 


1 + z 

10 


h i kpc. (6) 


A halo of given mass will collapse later in a universe 
with a non-zero cosmological constant, leading to the 
differences between the two sets of calculations presented 
here. The Einstein-de Sitter Universe calculations were 
first presented by BFPR, and here we have updated their 
diagram to reflect a modern cosmological model together 
with data for individual galaxies which, to the best of our 
knowledge, has not been done before. 

In order to compare with data, it is necessary to con¬ 
vert the M d yn defined by equation |T| to a total mass. 
A naive solution would be to assume that M to t ~ 


(3) 
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2M dyn (< R e ) but this ignores factors including variation 
in the mass to light ratio and the profile of the velocity 
dispersion. Instead, we use the ansatz that the total mass 
of a galaxy will be proportional to the dynamical mass 
within R e : 

M tot = OfMdyn (< R e ) (7) 

There will exist a scatter in a in the population of 
galaxies, and we may expect a to be a function of the 
total mass of the galaxy. However, if we approximate by 
considering a single value of a, we can then consider for 
what values of a the spherical collapse model discussed 
above is appropriate. Consider the example of a = 2, the 
naive value discussed above. The slope of the data is en¬ 
couraging, being in good agreement with the theoretical 
prediction. However, for a given mass, sigma is so large 
that the galaxies correspond to 4-sigma fluctuations in an 
EdS model and a 7-sigma fluctuation in ACDM. These 
are extremely unlikely fluctuations, not compatible with 
the observed number density of massive galaxies. 

In order to get a rough estimate for the level of fluctu¬ 
ations that correspond to massive galaxies, we use the 
morphologically segregated luminosity functions (LFs) 
from SDSS (Nakamura et al. 2003). Integrating the 
LF of E/SOs from L* to some bright upper limit, say 
20L* (the choice does not matter given the exponen¬ 
tial cutoff of the LF), we get a comoving number den¬ 
sity of ellipticals similar to those in our sample (see fig¬ 
ure The numbers from the SDSS LF give n(E/S0) = 
0.11 x 10 -2 /i 3 Mpc -3 . A similar integral for the LF of 
all galaxies in the sample in the range 0.005 <L/L* < 20 
gives the comoving number densities of galaxies, n(All) = 
9.02 x 10~ 2 /i 3 Mpc“ 3 . The ratio of the two number densi¬ 
ties (n(E/S0)/n(All) = 0.012) corresponds to a Gaussian 
fluctuation around 2.5-sigma. Hence, we can assume that 
a spherical collapse model from 2- to 3-sigma level fluc¬ 
tuations is compatible with the observed number density 
of early-type galaxies today. A choice of a = 20 produces 
a good correlation between the data and the predicted 
M-er relation for 3-sigma fluctuations given ACDM. 

It is interesting to note that this choice of a produces 
a value for the stellar mass to total mass ratio which 
is approximately half the cosmological baryon fraction; 
i.e. half the baryons in the halo are found in the form 
of stars. Models with large a are also consistent with 
the large mass to lig ht ratios found by st udies of massive 
ellipticals in x-rays llGriffiths et al.lll9 96'l. 

8. THE EVOLUTION OF THE DARK MATTER 

In order to further investigate the distributions of to¬ 
tal and stellar mass within the spherical collapse frame¬ 
work, we need to take into account the evolution of the 
dark matter as the baryons dissipate energy and sink 
towards the center of the halo. We use the procedure 
of adiabatic contraction of th e dark matter halo, follow¬ 
ing [Blume^^aOE^n (:1986b). We assume that the dark 
matter has an initial den sity profile of the form given by 
iNava.rro. Frenk <b White! : 1 997 1 (hereafter NFW). The 
concentration is defined as c = r v \ T /r s where r v - IT is the 
virial radius - usually considered to be the edge of the 
halo - and r s is the scale radius. This method is strictly 
valid only_jor halo particles moving on circular orbits. 
iGnedin et al.l (2004) show that the predicted density pro¬ 
file is more dense than that found in N-body simulations. 


Howe ver, the deviation is small iSellwood & McGauglj 
mi) and we use the method here as a first approxima¬ 
tion. 

The initial mass profile and the final dark matter and 
baryonic profiles are related by two equations: 

r [M g (< r) + M h (< r)] = rWb (< r»), (8) 

Mh (< r) = (1 — fcooi) Mi (< ri). (9) 

Mi (< r) is the initial mass enclosed within the radius,r. 
M g (< r) and Mh (< r) are the final baryonic and dark 
matter profiles respectively. f CO oi is the fraction of the 
system’s mass which is contained in baryons which cool 
to form the galaxy, and is related to a via the relation: 

f _ 2 M star {< R e ) 

U °° l ~ aM dyn (< R e ) [W) 

where the stellar mass within R, is assumed t o be half the 
cooled baryonic mass l|Padmanabhan et a 0 l200l l. For 
a = 20 and the ratio of M star (< R e ) to M dyn (< R e ) 
found in section 0] a value of f coo i ~ 0.08. This corre¬ 
sponds to half the cosmological baryon fraction, which 
suggests that half the halo’s baryons are in stars. We 
can also use this sort of argument to place a limit on a 
by requiring / coo 1 <1/6 (where 1/6 is the cosmological 
baryon density). This corresponds to a > 10. 

The baryons are assumed t o have a density profile of 
the form given by iHernauisfl (1199011 (which closely ap¬ 
proximates the de Vaucouleurs law, see §4). In this case 
the ad iabatic contra ction problem can be solved analyt¬ 
ically f Keetonl [2001 1. Each initial radius r* maps to a 
final radius r given by: 

fcooir 3 + (r + a) 2 [(1 - f coo i) r - r»] m, (< n) = 0, (11) 

where a = 0.551i? e is the scale radius of the Hernquist 
profile, and TOj(< n) = M t (< r;)/M tot . 

A halo of mass M tot in the monolithic picture consid¬ 
ered here will collapse at a redshift determined only by 
the cosmology and the amplitude of perturbations. The 
virial radius is then given by equation^] After choosing a 
value for a , / CO oi is given by equation^5]using the value of 
Af s tar (< 7? e )/Mdyn (< R e ) for each galaxy. For any con¬ 
centration, c, we then obtain the total mass inside the ini¬ 
tial radius (r*) that will evolve into a final radius r = R e . 
The adiabatic contraction prescription implies this mass 
is equal to the final mass observed within R e . Hence, 
we determine via a maximum likelihood estimator the 
concentration that satisfies M dyll (< R e ) = Mi (< n). 

8.1. Results 

Figure [7| shows the histogram of the concentration of 
the initial DM halos for our sample, for 2 and 3-sigma 
perturbations and a of 10 or 20. The figure overlays a 
Gaussian distribution corresponding to this result. The 
fit is remarkable: our simplistic model of spherical col¬ 
lapse and adiabatic contraction combined with a volume 
limited sample, results in a log-normal distribution of 
concentrations. 

This type of distributio n is also ob taine d in more de¬ 
tailed N-body simulations llBullock et all200lh . Further¬ 
more, the dispersion of our distribution is cr(logc) ~ 
0.2 — 0.25, in agreement with the simulations, which give 
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Concentration (initial) 


Fig. 7. — Concentration of the (NFW) dark matter halos before 
adiabatic contraction. The benchmark spherical collapse model, 
assuming LCDM and either 2- or 3-sigma perturbations is used. 
a = M to t/M dyn (< R e ) is either 20 or 10 as indicated. In all 
cases an approximately log-normal distribution of concentration is 
recovered. A higher order perturbation (3- instead or 2-sigma) or 
larger a reduces the mean value of the concentration. 


~ 0.18; notice that these simulations give values of con¬ 
centration higher than those we find even for 2-sigma 
fluctuations. 

Hierarchical models give concentrations around 4 only 
for objects which have recently undergone major merg¬ 
ers or for large (galaxy cluster mass) systems. The 
mean concentration for galaxies with mass M v ; r > 
1.5 x10 12 /i _1 Mq is typically 13.1 (Wechsler et al. 2002). 

As figure 0 shows, our model does reproduce large 
values of c for small values of a and 2-er fluctuations. 
However, a value a < 10 would shift the data points 
in figure 0 to the left, away from the predictions for 2- 
or 3-sigma fluctuations. We therefore conclude that the 
benchmark model predicts a mean concentration lower 
than that found in simulations. 

T hese results are_reminiscent of those of 

iCollister and Lahavl 1 )20051) . who studied the con¬ 
centration of galaxy clusters and groups in the 2dF 
galaxy redshift survey. They found concentrations 
from the galaxy data which were a factor of ~2 lower 
than the p redi c tions of dark matter simulations (e.g. 
iKravtsov et al.l (j2004h 1. albeit on larger scales than 
those considered here. Trujillo et al. (2006) present 
weak evidence for a change in concentration of massive 
elliptical galaxies (as derived from the size-mass relation) 
with redshift; further studies with larger samples should 
enable us to test their results with our methods. 

8.2. Galaxy formation and cooling 

In order to consider the baryonic component of the 
galaxies, it is convenient to plot baryon number density 
against temperature, following BFPR. Defining number 
density at present as: 


3M st ar(< Re) 
47 TRl/miH 


( 12 ) 



Fig. 8. — Galaxy formation a la BFPR Our sample galaxies are 
shown in a density-temperature diagram. The black dots represent 
the observed sample, whereas the grey dots “extrapolate” the sam¬ 
ple to the time of turnaround, before adiabatic contraction ensued. 
The grey and black lines are predictions for the collapse of uniform 
spheres corresponding to 2 and 3cr fluctuations, and are also shown 
before and after adiabatic contraction (AC). See text for details on 
the cooling curve; the dashed line corresponds to zero metallicity, 
the solid line to solar. 


and temperature as: 


T = nm H a 2 
kg 


(13) 


where fimg is the mean atomic weight (/r = 0.6 for a 
primordial mix of hydrogen and helium), and kg is the 
Boltzmann constant. The result is shown in figure[8] The 
sample is shown in two different ways: the black dots are 
the actual baryon densities (i.e. after adiabatic contrac¬ 
tion) using equation 1121 computed within the half-light 
radius. The grey dots represent baryon densities before 
adiabatic contraction, and are computed as the density 
corresponding to a baryon mass ( fcooi)M to t within a ra¬ 
dius R v i r . The grey and black lines are predictions for the 
spherical collapse model (ACDM), and are shown for 2 
and 3-sigma fluctuations before and after adiabatic con¬ 
traction. These models require a translation from total 
mass to baryon mass, for which we use the average value 
for fcool obtained from the sample, i.e. ( f coo i ) = 0.083. 
For the ‘before AC’ (grey) lines we use this mass and 
the virial radius definied by equation 0 Incidentally, we 
notice that the early-type galaxies indicated in the orig¬ 
inal cooling diagram in BFPR corresponds to our set of 
points before AC. As far as we know, this is the first 
time individual galaxies have been plotted in this space 
as previous work has relied on ‘schematic’ data. The 
(black) curves “after AC” also need to assume a ratio be¬ 
tween the virial radius and the initial radius (rf) which 
is mapped into R e . We use again the sample to find 
ri ~ 0.1R vir (see figure 0. The figure shows consistency 
with our previous results. 

The curves tracing td yn = t coo i uses the dynamical 
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r i/ r VIR 


Fig. 9.— Comparision of the virial radius r v ir and the initial ra¬ 
dius used in the adiabatic contraction procedure. For each galaxy 
Vi is chosen such that the mass within r{ before contraction cor¬ 
responds to the dynamical mass within R e after contraction. The 
top panel shows the correlation with the effective radius, R e . In 
larger galaxies, more of the halo is contracted into R e . 


timescale: 

tdyn = (14) 

where p is the total density of baryonic and dark matter 
and the cooling timescale is 

j. . 

tcool — . (lh) 

nKet 

where n is the number density of nuclei and A net is 
the net coo ling rate f or sola r (or zero) metallicity from 
ISutherland fc Dopital (1993). As expected, the galax¬ 
ies after contraction lie in the region of the plot where 
tcooi < tdyn- Although the approach is simplistic, assum¬ 
ing that galaxies are uniform spheres rather than using 
radial profiles, it is intriguing that 2-sigma fluctuations 
even for massive galaxies lie on the cooling curve for zero 
metallicity systems, and that adiabatic contraction pro¬ 
motes cooling of these systems. This result may be fol¬ 
lowed up with detailed simulations. 

9. CONCLUSIONS 

The sample of early-type galaxies selected by Bernardi 
et al. (2003a) from the Sloan Digitial Sky Survey has 
been used to investigate the properties of baryonic and 
dark matter within these systems. For the work pre¬ 
sented in this paper, the sample was further restricted to 
be volume limited, containing only massive elliptical sys¬ 
tems, leaving a total of more than 2000 systems. Using 
photometry from Sloan Data Release 4 and the synthetic 
stellar populations of Bruzual & Chariot (2003) we cal¬ 
culated the mass in stars of each system. 

We recover the well-known Faber-Jackson relation be¬ 
tween luminosity and observed velocity dispersion and 


also find a tight correlation between stellar mass and 
velocity dispersion, which we have called the ‘bary¬ 
onic Faber-Jackson relation’. We find a slope of 
A/star (< Re) oc g 2 ' 33=1=0 ' 1 8 (in approp ri ate u nits), less 
steep than previous work (iThomas et alJl2005l) based on 
a study of an order of magnitude fewer galaxies. 

In addition to determining the stellar mass, it is possi¬ 
ble to determine a dynamical mass via the velocity dis¬ 
persion and the radius of the system. Comparing the 
two estimates of mass reveals that dark matter is under- 
abundant in the centre of such galaxies when compared 
to the cosmological baryon fraction. However, our results 
are in good agreement with observers who have claimed 
that the inner regions of massive elipticals are baryon 
dominated. 

In order to examine whether this large data set is 
sufficient to discriminate between galaxy formation via 
monolithic collapse or the merging of small dark mat¬ 
ter haloes, the data were compared to a simple model of 
spherical collapse. It was necessary to assume that the 
total mass of the systems under study is proportional to 
the mass within the baryon-rich core (within R e ), but 
given such an assumption we find that for a ~ 20 (where 
a = M t ot/Md y n (< R e )), the data were found to be con¬ 
sistent with simple monolithic collapse of 3-sigma petur- 
bations (assuming ACDM cosmology). 

Interestingly, this choice of a also produces a value for 
the fraction of the mass which lies in the form of cooled 
baryons which is consistent with other work. A choice of 
20 seems large, but corresponds to the stellar mass con¬ 
tributing half of the cosmological baryon fraction. Given 
the other reservoirs of baryons in these systems, this is a 
reasonable value. 

Having motivated our choice of a, we used t h e adia¬ 
batic co ntraction method first due to iBlumenthal et all 
( 1986b ) to calculate the response of the dark matter to 
the presence of the observed baryons. We assume that 
the dark matter has an initial NFW profile, and produce 
the distribution of initial concentration shown in figure 
0 The mean concentration is lower than that found in 
dissipation-less N-body simulations. 

In summary, we have shown that the B03 sample of 
early-type galaxies is a powerful resource for studying the 
distribution of matter within these objects. However, de¬ 
spite this excellent data set we cannot rule out spherical 
collapse as a mode of formation for massive ellipticals. 
We find that the central regions of the galaxies, which 
are baryon dominated, must represent only a small frac¬ 
tion of the total mass of the system. In future work, we 
will extend our analysis to a new large sample of early- 
type systems, and undertake a detailed comparision of 
simulation results in order to understand the differences 
between observations and semi-analytic models. 
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